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ABSTRACT 

We present integral field spectroscopy of two selected zones in the Orion Nebula ob- 
tained with the Potsdam Multi- Aperture Spectrophotometer (PMAS), covering the 
optical spectral range from 3500 to 7200 A and with a spatial resolution of 1" . The 
observed zones are located on the prominent Bright Bar and on the brightest area at 
the northeast of the Orion South cloud, both containing remarkable ionization fronts. 
We obtain maps of emission line fluxes and ratios, electron density and temperatures, 
and chemical abundances. We study the ionization structure and morphology of both 
fields, which ionization fronts show different inclination angles with respect to the 
plane of the sky. We find that the maps of electron density, 0+/H+ and 0/H ratios 
show a rather similar structure. We interpret this as produced by the strong depen- 
dence on density of the [O ll] lines used to derive the 0+ abundance, and that our 
nominal values of electron density -derived from the [S ll] line ratio- may be slightly 
higher than the appropriate value for the 0+ zone. We measure the faint recombina- 
tion lines of O ll in the field at the northeast of the Orion South cloud allowing us 
to explore the so-called abundance discrepancy problem. We find a rather constant 
abundance discrepancy across the field and a mean value similar to that determined 
in other areas of the Orion Nebula, indicating that the particular physical conditions 
of this ionization front do not contribute to this discrepancy. 
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1 INTRODUCTION 

The analysis of H II regions is currently our best tool to study 
the physical conditions and chemical abundances of the in- 
terstellar medium. Because of its proximity and high surface 
brightness, the Orion Nebula (M42) is one of the best stud- 
ied H II regions of our galaxy. Although the Orion Nebula 
has a large extension for more than one half degree, most of 
the radiation comes from the inner part, the so-called Huy- 
gens region. This is an active star-formation region ionized 
by a group of four massive stars known as the Trapezium 
cluster, where O^Ori C (07 V; Simon-Diaz et al. 2006) is the 
main ionizing source and, therefore, responsible of the bulk 

* Based on observations collected at the Centro Astronomico His- 
pano Aleman (CAHA) at Calar Alto, operated jointly by the 
Max-Planck Institut fiir Astronomic and the Instituto de As- 
trofisica de Andalucia (CSIC). 
t E-mail: amd@ifa.hawaii.edu/adalfis@gmail.com 



of nebular emission from the main ionization front (MIF). 
Physical, chemical, kinematical and structural properties for 
this H II region have been studied in depth using different ob- 
servational techniques {e.g. Hester et al. 1991; Pogge et al. 
1992; Walter et al. 1992; Esteban et al. 1998; Bally et al. 
2000; O'Dell 2001; O'Dell et al. 2003; Doi et al. 2004; Es- 
teban et al. 2004; Garcia-Diaz & Henney 2007; Garcia-Diaz 
et al. 2008; O'Dell & Henney 2008; O'Dell et al. 2009) as 
well as photoionization models {e.g. Zuckerman 1973; Bal- 
ick et al. 1974; Rubin et al. 1991; Baldwin et al. 1991; Wen 
& O'Dell 1995; O'Dell 2001). The presence of many different 
morphological features such as protoplanetary discs, Herbig- 
Haro objects or bright bars makes the Orion Nebula the best 
laboratory for analysing photoionization phenomena. 

They study of ionization fronts is one of these phenom- 
ena of special interest because they are transitions between 
the fully ionized gas and its neutral surroundings, where we 
expect density and temperature gradients. Probably, there 
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Figure 1. HST image of the central part of the Orion Nebula, which combines WFPC2 images taken in different narrow-band filters 
(O'Dell & Wong 1996). The white squares correspond to the fields of view of PMAS IFU used, covering the Bright Bar and the northeast 
edge of the Orion-S cloud (NE-Orion-S). The black and white images show enlargements of the Ho; emission in our PMAS fields, rebinned 
to 160x160 pixels and smoothed using a Gaussian filter. 



is an ionization front in almost every direction towards the 
Orion Nebula, though only certain geometries allow us to 
explore the ionization structure of these features. The clear- 
est example in the Orion Nebula is the well-known Bright 
Bar (hereinafter, BB). It is one of the most conspicuous fea- 
tures of the Orion Nebula, which is seen as an elongated 
structure in both ionized gas and molecular emission. It is 
located at about 111'' to the southeast of 6'^Ori C. The BB 
appears specially bright in emission lines of low ionization 
species such as [N ii] or [S ii] [e.g. Hester 1991; Pogge et al. 
1992; Walter et al. 1993; Garcia-Diaz & Henney 2007). On 
the basis of this fact, Balick et al. (1974) proposed that the 
BB is an escarpment in the MIF. Multi- wavelength studies 
from infrared to radio spectral domains as well as theoretical 
modeling {e.g. Tielens et al. 1993; Tauber et al. 1994; van 
der Werf et al. 1996; Marconi et al. 1998; Young Owl et al. 
2000; Kassis et al. 2006; van der Wiel et al. 2009; Pellegrini 
et al. 2009, and references therein) support this proposal, 
where the MIF changes from a face-on to an edge-on geom- 
etry. Pellegrini et al. (2009) modeled the emission spectra 
along the transition from to H° and H2 regions. They 
adopt a thin slab geometry for the BB, which extends 0.115 
pc and is inclined 7° with respect to the line of sight, finding 
a good agreement between the computed optical spectrum 
and that observed by Baldwin et al. (1991). 

Besides the BB, there are other remarkable areas in 
the Orion Nebula. Recently, Garcia-Diaz & Henney (2007) 



have identified several examples of compact bars toward the 
core of the Huygens region, in addition to the brighter bars 
previously reported in the literature (O'Dell & Yusef-Zadeh 
2000). This is the case of the brightest zone of the Huygens 
region (see Fig. 1), located at the northeast at the north- 
east edge of the Orion South (Orion-S) cloud (hereinafter, 
NE-Orion-S). This region is located about 30'' to the south- 
west of ^^Ori C. In several observational studies {e.g. Pogge 
et al. 1992; Walter et al. 1993; Garcia-Diaz & Henney 2007) 
the NE-Orion-S appears as an extended structure, which is 
very bright in [N ii] and [S ii] emission lines as in the case 
of the BB. Recently, O'Dell et al. (2009) have studied in- 
depth the three dimensional structure of the whole Orion-S 
concluding that it is a cloud ionized from one side facing to 
^^Ori C and suspended within the main body of the Orion 
Nebula in front of the MIF. Therefore, this cloud as well as 
its NE-Orion-S edge are ionization fronts with a given un- 
known inclination, though it should not be as tilted as the 
BB (O'Dell et al. 2009). 

As we will see in the next paragraphs, there have been 
a number of works devoted to study the spatial distribu- 
tions of physical conditions and chemical abundances in the 
Orion Nebula, but this is the first one specifically dedicated 
to these ionization fronts -BB and NE-Orion-S- making use 
of integral field spectroscopy in the optical spectral range at 
spatial scales of about l". Walter et al. (1992) determined ra- 
dial gradients of the physical conditions from the analysis of 
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Table 1. Journal of observations. 
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" Coordinates of the field centre (J2000.0). 



22 different areas of the nebula obtained from long-slit spec- 
trophotometry. Pogge et al. (1992) and Walter et al. (1993) 
studied the ionization structure of the entire Huygens region 
by means of Fabry-Perot and CCD imaging, respectively. 
Pogge et al. (1992) obtained the first electron density map, 
finding the maximum values precisely at the NE-Orion-S 
position. More recently, O'Dell et al. (2003) carried out a 
high-spatial resolution map of the electron temperature (de- 
rived from the [O III] line ratio) centered at the southwest of 
the Trapezium cluster from fiux-calibrated narrow-band im- 
ages taken with the Hubble Space Telescope (HST). Rubin 
et al. (2003) also analysed long-slit data taken with STIS 
at the HST, obtaining the spatial distributions of electron 
temperatures (derived from the [O III] and [N ll] line ratios) 
along 4 slit positions. Using long-slit spectroscopy at spa- 
tial scales of l'!2 and covering the Huygens region with 5 
slit positions, Mesa-Delgado et al. (2008) also obtained the 
spatial distributions of emission lines fluxes, physical con- 
ditions and chemical abundances, specially of those derived 
from faint recombination lines of heavy-element ions. There 
are just a few works using integral field spectroscopy in the 
Orion Nebula. Sanchez et al. (2007) obtained a mosaic of the 
Orion Nebula from integral field spectroscopy with a spatial 
resolution of 2'.' 7 but with very short exposure times. The 
electron temperature map obtained by these authors from 
the [N ll] line ratio shows spatial variations and the elec- 
tron density map is very rich in substructures. There are 
only three additional studies available in the literature that 
used this observational technique in small fields of the Orion 
Nebula, though they are focused in the analysis of certain 
morphological structures. These are the cases of Vasconcelos 
et al. (2005) and Tsamis et al. (2011), who study the proto- 
planetary disc LV2, and Mesa-Delgado et al. (2009), focused 
in the study of the prominent Herbig-Haro object HH 202. 

In §2 we describe the observations made with PMAS 
and the reduction procedure. In §3 we describe the emission 
line measurements and the reddening correction of the spec- 
tra. In §4 we explain the calculations of physical conditions 
and chemical abundances. In §5 we present the spatial dis- 
tributions of the flux of several emissions lines and emission 
line ratios, physical conditions and chemical abundances. In 
§6 we discuss the structure of the two ionization fronts stud- 
ied as well as the so-called abundance discrepancy problem. 
Finally, in §7 we summarize our main conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

In Fig. 1 we show the two fields on the Orion Nebula ob- 
served in service time on 2008 December 19-20 at Calar Alto 
Observatory (Almeria, Spain). Integral field spectroscopy 
was performed in these fields using the Potsdam Multi- 
Aperture Spectrophotometer (PMAS, Roth et al. 2005) at 



the 3.5-m Telescope and its standard lens array integral field 
unit (IFU) of 16''xl6'' field of view (FoV) with a spaxel size 
of l''xl''. The optical range from 3500 to 7200 A was cov- 
ered using the V600 grating and two rotator angles (—72° 
for the blue spectra and —68° for the red one). This grating 
allowed us to achieve an effective spectral resolution of about 
3.6 A. The total integration times for the blue and red spec- 
tra as well as the central coordinates of each field are pre- 
sented in Table 1. Additional short exposures of 10 seconds 
were also taken in order to avoid saturation of the brightest 
emission lines. Calibration images were also obtained dur- 
ing the night: HgNe arc lamps for the wavelength calibration 
and continuum lamps needed to extract the 256 individual 
spectra on the CCD. The PMAS data were reduced (bias 
subtraction, spectra extractions and flat-fielding correction) 
as well as wavelength and flux calibrated using the IRAF^ re- 
duction package specked (see Mesa-Delgado et al. 2009, for 
more details on the reduction process). The data were cal- 
ibrated using the spectrophotometric standard stars Feige 
100, Feige 34, Feige 110 and G 191-B2B (Oke 1990) and 
the error of the absolute flux calibration was about 5%. The 
night reports indicated the presence of thin clouds during 
the observations of the BB, and this could be the most likely 
reason of a disagreement of about 15% in the absolute flux 
measurements between the short and long time exposures of 
the BB in the red range. We corrected for this effect and we 
do not expect that the final results are substantially affected 
because we only use emission line ratios in our analyses. This 
correction was also considered in the calculation of the flux 
measurement errors (see §3.1) for the spectra of this range 
and field. In the blue range of the BB field we found a dif- 
ference between the short and long exposure spectra lower 
than 3% (even lower than the nominal error in the flux cal- 
ibration). For the NE-Orion-S, the differences between the 
short and large exposures were also lower than 3% along 
the whole wavelength coverage in both blue and red spec- 
tral ranges. The seeing had an average value of about l'.'5 
during the observations. 

The integral fleld data cubes are subject to the effects 
of the differential atmospheric refraction (DAR) along the 
direction of the parallactic angle. This effect was corrected 
using our own IDL routines (similar to the previous ones writ- 
ten by P.T. Wallace for STARLINK and Walsh k Roy 1990) 
implementing the method outlined by Hohenkerk & Sinclair 
(1985) (see also §3.281 of Seidelmann 1992) The procedure 
calculates the fractional spaxel shifts for each monochro- 
matic image with respect to a given wavelength. The max- 
imum DAR shift was obtained for the BB field reaching 
values of Aa = 0'.'62 and = 1'.'16 between [O ii] 3727 

^ IRAF is distributed by NO AO, which is operated by AURA, 
under cooperative agreement with NSF. 
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A and The reference wavelengths were H/3 and Ka for 
the blue and red data cube, respectively. These emission 
lines were also used to align the blue and red data cubes 
calculating the right spaxel offsets that produce the maxi- 
mum correlation between both H/3 and Ha emission maps. 
For the two fields, these offsets reached values lower than 
0'/3 in a and S. Finally, the data cubes are corrected using 
the BILINEAR function of IDL, the calculated fractional shifts 
and the alignment offsets. As result of all this process, the 
maximum coincident FoV resulting in the whole wavelength 
range is reduced to lA" xlA" . We tested the DAR correction 
observing the alignment of localized features at different H I 
lines finding good agreement in all cases. 
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Figure 2. Section of a PMAS spectrum around O ll lines of mul- 
tiplet 1 corresponding to the spaxel position ( — 1:— 2,5:6) of the 
NE-Orion-S field. The O ll feature at about 4650 A corresponds 
to the blend of O ii 4649, 4651 A lines. The O n+N in+N ii fea- 
ture corresponds to the blend of O ii 4639, 4642 A, N ill 4641, 
4642 A, and N II 4643 A lines. 



3 LINE MEASUREMENTS AND REDDENING 
CORRECTION 

3.1 Line intensity measurements 

We measured the intensity of all the emission lines needed 
to perform our analysis: hydrogen B aimer lines (from Ha 
to H12), which are used to compute the reddening correc- 
tion and verify the DAR correction; coUisionally excited lines 
(CELs) of different species (see §4.2), which are needed to 
compute physical conditions and chemical abundances; and, 
finally, the faint C II and O II lines, which are used to de- 
rive the ionic abundances from recombination lines (RLs). 
It should be mentioned that both RLs were only detected 
in the NE-Orion-S field. These faint lines were not detected 
in the BB due to its lower surface brightness and lower ion- 
ization degree. In Fig. 2, we present the spectra of a repre- 
sentative spaxel, ( — 1:— 2,5:6)^, of the NE-Orion-S, where we 
can see the blend of O II 4649, 4651 A lines used to derive 
the 0^+/H+ ratio. 

The measurement of the emission line fluxes were per- 
formed applying a single -or a multiple, in the cases where 
it was necessary- Gaussian profile fit procedure between two 
given limits and over the local continuum making use of the 
SPLOT routine of IRAF. All measurements of the bright lines 
were made with our own scripts to automatize the process. 
However, the fiux of the faint C II 4267 A and O ii 4650 
A lines of the NE-Orion-S field were measured by hand for 
each individual spaxel. We have determined the errors in the 
flux measurements following the criteria deflned by Mesa- 
Delgado et al. (2008). The flnal error for each emission line 
was the result of the quadratic sum of the flux measurement 
error and the flux calibration error. As in our previous pub- 
lications, we defined three criteria to avoid spurious weak 
line measurements in the automatic procedure and to dis- 
criminate between real features and noise: 1) line intensity 
peak over 2.5 times the sigma of the local continuum; 2) 
FWHM(H i)/1.5 < FWHM(A) < 1.5xFWHM(H i); and 3) 
F(A) > 0.0001 X F(H^). 



3.2 Reddening correction 

The reddening coefficient, c(H/3), was obtained by fitting 
the observed H7/H/3 and H5/H^ ratios to the theoretical 
ones predicted by Storey & Hummer (1995) for Ue = 1000 
cm~^ and Te = 10000 K. In a previous step, we considered 
more H l line ratios in the calculations, but they only added 
more dispersion to the final c(H/3) value because the rest 
of available lines were either rather faint (Hll and HI 2), 
or blended with other bright lines (the case of H7, H8 and 
H9). We used the updated reddening function, /(A), normal- 
ized to H/3 for the Orion Nebula [Ry — 5.5; Blagrave et al. 
2007) . The use of this extinction law instead of the classical 
one of Costero & Peimbert (1970) produces c(H^) values 
about 0.1 dex higher and dereddened fiuxes with respect to 
Hy5 about 3% lower for lines in the range 5000 to 7500 A, 
4% higher for wavelengths below than 5000 A and 2% for 
wavelengths above 7500 A. The final adopted c(H/3) value 
for each spaxel is an average of the individual values derived 
from each B aimer line ratio weighted by their corresponding 
uncertainties. The typical error of c(H/3) is about 0.20 dex 
for each spaxel. All emission line fiuxes of a given spectrum 
were normalized to H/3 and Ha for the blue and red range, 
respectively. To produce a final homogeneous set of line flux 
ratios, all of them were re-scaled to H/3 after the reddening 
correction. The re-scaling factor used in the red spectra was 
the theoretical Yia/Yifi ratio for the physical conditions of 
Te= 10000 K and ne= 1000 cm"^ 

The resulting extinction maps are shown in Figs 3a and 
3b for the BB and NE-Orion-S, respectively. The extinc- 
tion coefficients vary approximately from 0.4 to 1.0 dex for 
the two fields. In the case of the BB, the highest values, 
0.9 and 1.0 dex, are located at the southeast corner of the 
field, while the lowest values are found at the position of the 
ionization front in the middle of the field. The NE-Orion-S 

Hereinafter, the localization of an individual spaxel or multiple 
spaxels will be given as (X1:X2,Y1:Y2) according to the coordi- 
nate system used in the maps presented in next sections. 
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Figure 3. Spatial distribution maps of the extinction coefficient, c(H/3), with Ho; contours overplotted for the two PMAS fields: a) Bright 
Bar and b) NE-Orion-S. 



presents values between 0.5 to 0.7 dex in most of the field 
and somewhat larger values -between 0.7 and 0.9 dex- at 
the east edge of the observed area. We have estimated the 
mean c(H/3) value of each field as well as the standard de- 
viation, which amounts to 0.7=b0.2 and O.GzbO.l dex for the 
BE and NE-Orion-S, respectively. These mean values are in 
agreement with other determinations available in the liter- 
ature. O'Deh & Yusef-Zadeh (2000) obtained c(H/3) maps 
and derived consistent c(H/3) values both from the 
line ratio by using calibrated HST images and from radio 
to optical surface brightness ratio. Using their extinction 
maps, we have found c(H/3) values ranging from 0.4 to 0.6 
for the BB and 0.4 to 0.7 for the NE-Orion-S at the precise 
positions of our fields. These extinction values of O'Dell & 
Yusef-Zadeh (2000) were determined using the classical ex- 
tinction law of Costero & Peimbert (1970) and, as we com- 
mented before, should be increased by about 0.1 dex to be 
compared with our determinations. 



4 DETERMINATIONS OF PHYSICAL 
CONDITIONS AND CHEMICAL 
ABUNDANCES 

4.1 Physical conditions 

We have used the temden task of the nebular package 
of IRAF (Shaw & Dufour 1995) with updated atomic data 
(see Garcia- Rojas et al. 2009), to determine the physical 
conditions, electron densities, Ue, and temperatures, Te, for 
each spaxel of our PMAS fields from the usual GEL ratios 
-[S II] 6717/6731 A for n^, and [O iii] (4959+5007) /4363 
A and [N ii] (6548+6584) /5755 A for Te. For these calcu- 
lations, we have followed the same procedure described in 
Mesa-Delgado et al. (2008) assuming an initial Te= 10000 
K to obtain a first approximation of Ue, then we calculate 
Te([0 III]) and Te([N ii]), and iterate until convergence. We 
have not considered the contribution by recombination in 
the observed flux of the auroral line [N ii] 5755 A to derive 
Te([N ii]), because it is very small for the conditions of the 
Orion Nebula (e.g. Esteban et al. 2004). The typical uncer- 
tainties in the density maps of the BB range from 1000 to 
1200 cm~^. For the NE-Orion-S field, the typical error in the 
density distribution amounts to 4000 cm~^. This higher er- 
ror is due to the much higher densities at this particular area 
of the nebula {e.g. Pogge et al. 1992; Wen & O'Deh 1995; 



Mesa-Delgado et al. 2008), that makes the [S ii] indicator to 
be about the high density limit. The typical error for elec- 
tron temperatures range from 500 to 1000 K for Te([N ii]) 
and from 100 to 200 K in the case of Te([0 ill]). 

4.2 Chemical abundances 

We have derived abundances for several ions from GELs: 
N+, 0+, 0^+, S+, S^+, Ne2+ and At^+ making use of the 
IRAF task IONIC of the package nebular. We have assumed a 
two-zone scheme and adopting Te([N ii]) for low ionization 
potential ions (N+, 0+ and S+) and Te([0 ill]) for high 
ionization potential ones (O^^, S^^, Ne^^ and Ar^^). We 
have computed the ionic abundance errors as the quadratic 
sum of independent contributions from Ue , Te , and line flux 
uncertainties. For ionic and total abundances of oxygen, the 
average uncertainties are 0.12-0.20 dex for 0^/H+, 0.03 dex 
for 0^+/H+ and 0.10-0.15 dex for 0/H. 

In addition, we have also detected and measured pure 
RLs of O II and G ii in the majority of spaxels covering 
the NE-Orion-S field (see Fig. 2). RLs have the advantage 
that their relative fluxes with respect to H i lines depend 
weakly on Te and Ue, avoiding the problem of the possible 
presence of temperature and/or density fluctuations. These 
fluctuations can be actually affecting the abundance deter- 
minations from GELs, whose emissivities strongly depend on 
the physical conditions of the ionized gas (Peimbert 1967). 

The spectral resolution of our spectra does not allow to 
separate the blend of the two brightest individual lines of 
multiplet 1 of O II at 4649 and 4651 A. In this situation, the 
combination of equations (1) and (4) by Peimbert & Peim- 
bert (2005) should be used to determine the O^^ abundance 
from these RLs. These expressions permit to correct the line 
fluxes from NLTE effects in the relative intensity of the indi- 
vidual lines of the multiplet, though they are rather small in 
the Orion Nebula due to its relatively high density. There- 
fore, we have calculated the 0^^/H^ ratio from RLs as: 

O^^Am aeff{Hl3) /(Ml O II) 

H+ 4861 ^ aeff (Ml) ^ I{H^) ' ^ ^ 

where aeff{Hl3) and aeff (Ml) are the effective recombina- 
tion coefficients for H/3 and for the O ii multiplet 1, respec- 
tively, and Ami = 4651.5 A the representative mean wave- 
length of the whole multiplet. We have also calculated the 
G^+ /H+ ratio from the flux of G ii 4267 A RL and using 
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Figure 4. Spatial distributions of emission line fluxes in units of erg cm ^ s -"^ of the main ionization stages of oxygen at the conditions 
of the Orion Nebula: [O l] 6300 A (top), [O ll] 3727 A (middle) and [O ill] 4959 A (bottom). The maps correspond to the Bright Bar 
(left) and NE-Orion-S (right) flelds. Ho; contours are overplotted in all maps. 



an analogous equation to (1) particularized for the case of 
this ion. According to the two-zone scheme, we have adopted 
Te([0 III]) to calculate the abundances of both ions, O^^ and 
C^^. We have used the effective recombination coefficients 
of Storey 1994 to derive the O^^ abundances (assuming LS 
coupling) and Davey et al. 2000 for C^^ abundances. The 
typical errors in the abundances from RLs amount to 0.08- 
0.10 dex for C^+/H+ and 0.10-0.12 dex for 0^+/H+. 



5 RESULTS: SPATIAL DISTRIBUTION MAPS 

In this section we analyse the spatial distributions of several 
emission line fluxes, emission line ratios, electron densities 
and temperatures as well as abundances for several repre- 
sentative ions for each PMAS field. 



5.1 Emission line fluxes and excitation ratios 

In Fig. 4, we present the spatial distribution maps of [O i] 
6300 A, [O ii] 3727 A and [O iii] 4959 A emission hne fluxes 
for both the BB (Figs 4a-4c) and NE-Orion-S (Figs 4d- 
4f). These maps show the distribution of the different ionic 
species of oxygen that can be found in the Orion Nebula: 0°, 
0+ and 0^+. As it can be seen, the ionization stratification 
of the two ionization fronts is clearly resolved in the maps. 
However, these maps are really showing a combination of 



the ionization stratification and the density distribution of 
the gas because the emissivity of a line is proportional -at 
first order- to Ue^. It is very illustrative to compare these 
maps with those of the [O i]/Ha, [N ii]/Ha and [O iii]/H/3 
line ratios shown in Fig. 5. Line ratios are directly related to 
the local ionization degree of the gas and do not depend in 
a direct manner on density. In Fig. 5 we represent [N ii]/Ha 
instead of [O ii]/H/3 because the [N ii]/Ha ratio is the most 
commonly used indicator for the analysis of ionization fronts 
and shows a very similar spatial distribution to [O ii]/H/3. 
As it is explained in O'Dell et al. (2009), the [N ii]/Ha ratio 
is a good indicator to detect when the MIF is tilted and, 
therefore, this is the reason that the BB is so well defined 
in Fig. 5. 

As it has been said before. Figs 4-5 show the ionization 
stratification at and around the ionization fronts . In both 
figures, it is evident that the areas with the higher [O l] 
6300 A line emission and [O i]/Ha ratio are narrower than 
those of the other lines represented in the figures, and their 
maxima tend to be also located farther away with respect to 
^^Ori C (see Fig. 1). We can notice some differences between 
the spatial distributions of line fluxes and ratios for the BB 
and NE-Orion-S. Firstly, the widths of the zones occupied 
by the tracers of the ionization fronts, [O i] and [N ii], are 
narrower for the BB than for the NE-Orion-S. This indicates 
that the inclination angles of the planes of the fronts with 
respect to the line of sight of both features are different. In 
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Figure 5. Spatial distribution maps of selected emission line ratios for the Bright Bar (left) and NE-Orion-S (right). The sum of [N ll] 
6548, 6584 A and [O iii] 4959, 5007 A nebular hues have been used to derive the [N ii]/Hck; (middle) and [O m]/¥L/3 (bottom) ratios, 
respectively. Ho; contours are overplotted in all maps. 



fact, as we commented in §1, this was already suggested by 
O'Dell et al. (2009). The BB is located at a border of the 
Huygens region where the MIF is an edge-on escarpment, 
where its plane is tilted about 7° with respectto our line of 
sight (Pellegrini et al. 2009). The situation is not so clear 
in the case of the NE-Orion-S. It is a portion the Orion- 
S cloud, which is suspended in front of the MIF with an 
unknown inclination and can be interpreted as a bump in the 
otherwise concave surface of the MIF at the Huygens region 
(O'Dell et al. 2009). Secondly, another difference between 
the BB and NE-Orion-S can be seen comparing the maps 
shown in Figs 4 and 5. In the case of the BB, Fig. 4b) shows 
that the maxima of [O ii] 3727 A and Ha emission coincide. 
However, in the case of the NE-Orion-S, Fig. 4f) indicates 
that the maximum of Ha emission coincides with [O ill] 
4959 A. This difference is due to the NE-Orion-S area is 
nearer to ^^Ori C than the BB and that while few photons 
capable of producing O^^ are reaching the BB, NE-Orion-S 
is illuminated by a significant fraction of such high-energy 
photons. This is clearly confirmed in Fig. 5c), where the 
maximum of the [O lll]/H/3 line ratio is well outside the BB 
and towards the ionizing source of the nebula. Interestingly, 
Fig. 5f) indicates that in the case of the NE-Orion-S, the 
zone of the highest [O iii]/H/3 line ratio is just at the north 
edge of the field. This area coincides with a filament, which 
seems to be out of Orion-S cloud itself. Its higher excitation 



suggests that this filament should be nearer to 6'^Ori C and 
perhaps belongs to the MIF. 



5.2 Physical conditions 

In Fig. 6 we present the physical conditions for the two ob- 
served fields: ne([S ii]), Te([N ii]), Te([0 ill]) and the tem- 
perature ratio Te([N Ii])/Te([0 ill]). 

Both fields show Ue spatial variations of the order of a 
factor 2. The densities in the BB range from 3000 to 5500 
cm~^, but larger values are found in the NE-Orion-S, from 
8000 to 16000 cm~^. In the BB field, maximum densities are 
found within the Ha contours, while in the NE-Orion-S these 
are displaced with respect to the maximum of H i emission. 
The average densities of the ionization fronts are 4370=b580 
cm"^ for the BB and 10900±1900 cm"^ for the NE-Orion-S. 
In the case of the NE-Orion-S field, the highest density val- 
ues are about the high density limit of the [S ii] line ratio in- 
dicator and, therefore, they may not correspond to the true 
densities. Fortunately, we have detected at least 5 [Fe ill] 
lines of the 3F multiplet ([Fe ill] 4607, 4658, 4702, 4734 and 
4755 A) and one line of the 2F multiplet ([Fe ill] 4881 A) at 
those spaxels of the NE-Orion-S field showing the ne([S ii]) 
maxima. These lines permit to get an additional determina- 
tion of the density, ne([Fe ill]), which range of validity has 
an upper limit above that of the [S ii] indicator. To carry 
out these calculations we have used a 34-level atom using 
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Figure 6. Spatial distributions of physical conditions for the Bright Bar (left column) and NE-Orion-S (right column). For all fields, we 
present ne([S ll]) (top row) in units of cm~^, Te([N ll]) (second row from top) and Te([0 ill]) (third row from top) in units of K, and 
the ratio of both temperatures (bottom row). Ho; contours are overplotted in all maps. 



the collision strengths from Zhang (1996) and the transition 
probabilities of Quinet (1996) and Johansson et al. (2000). 
The calculations of ne([Fe ill]) give densities varying from 
10000 to 18000 cm~^, even though they have large uncer- 
tainties (between 3000 and 6000 cm~^) due to the faintness 
of the [Fe ill] lines. The differences among the ne([Fe ill]) and 
ne([S ll]) determined for a given spaxel amount from 1% to 
30%, and it is clear that these differences are within the un- 
certainties of both indicators. Therefore, the high ne([S ii]) 
we obtain in Fig. 6e) should be considered at least close to 
the true electron density of the ionized gas at the NE-Orion- 
S (further discussion of these results can be found in §5.3). 

The spatial distributions of Te([N ii]) show a similar 
qualitative behaviour in both ionization fronts (Figs. 6b and 
6f). The pattern seems to be almost inverse to the density 
distributions reaching the minimum temperatures approxi- 
mately at the maximum densities. In the case of the BB, 
the range of variation of Te([N ii]) along the field is rather 
narrow reaching the maximum temperatures -about 9500 
K- just at the outer edge of the bar. In the case of the NE- 
Orion-S field, the range of values of Te([N ii]) is wider and 
the temperatures are always higher than in the BB. The 



maximum values of temperature are also found at the outer 
edge of the ionization fronts. 

The spatial distributions of Te([0 ill]) (Figs 6c and 6g) 
are rather featureless with mean values of about 8620 and 
8670 K and variations of the order 600 and 400 K in the 
BB and NE-Orion-S fields, respectively. In both fields, the 
maximum values of Te([0 ill]) tend to be located within the 
H I emission. In the NE-Orion-S, we can see a zone of higher 
values located at the southwest corner of the field, though 
these values have also higher uncertainties -about 200-250 
K. 

In Figs 6d) and 6h) we present the Te([N Ii])/Te([0 ill]) 
ratio for both fields. As we can see, Te([N ii]) is always 
slightly higher than Te([0 ill]). This is the expected be- 
haviour in an ionization bounded nebula as a consequence 
of the hardening of the incident radiation field due to pho- 
toelectric absorption in the low ionization zones where the 
[N ii] lines are emitted. This ratio tends to be higher in the 
outer zones of both fronts. 

We have compared our determinations of the physical 
conditions with previous ones available in the literature for 
the same or nearby areas. These are the cases of: the long- 
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Figure 7. Spatial distributions of the 0"^/H+, 0^"^/H"^ and O/H ratios -in the usual logarithmic scale- for the Bright Bar (left column) 
and NE-Orion-S (right column). Ho; contours are overplotted in all maps. 
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Figure 8. Spatial distributions of and 0"^+ abundances determined from recombination lines for the NE-Orion-S. Ha contours are 
overplotted. 



slit positions 1, 3, 4 and 5 analysed by Mesa-Delgado et al. 
(2008), which cover a small area of the Orion-S cloud and the 
BE; the long-slit positions 2 and 4 studied by Rubin et al. 
(2003), which cover some parts of both ionization fronts; sev- 
eral of the slit positions observed by Baldwin et al. (1991); 
the results presented by Pogge et al. (1992) for the whole 
Huygens region; and the results from the model of Wen & 
O'Dell (1995). In general, we have found a good agreement 
between the results obtained in these papers and our deter- 
minations, reproducing some of the patterns discussed in the 
previous paragraphs. However, we have found a slight dis- 
agreement with the Ue determinations of Pogge et al. (1992) 



for the BE. These authors analysed a spatial profile of Ue 
crossing the EE in an area located about 40^^ to the south- 
west of our field and with a position angle of 150°. Their pro- 
file reaches densities of about 3200 cm^'^-lOOO cm~^ lower 
than our determination- but in agreement with the results 
by Garcia-Diaz & Henney (2007) along a similar slit po- 
sition. On the other hand, our density determination is in 
agreement with the model results of Wen & O'Dell (1995). 
In their figure 4a, they present a spatial profile of the den- 
sity with a position angle of 322°, which precisely crosses 
our EE field, finding densities of about 5000 cm~^, similar 
to our values. Therefore, the small differences found here be- 
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tween our physical conditions and those reported by other 
authors seem to be related to the presence of local varia- 
tions of the physical conditions in the BB. Wen & O'Dell 
(1995), in their figure 5a, also presented a spatial profile of 
the density with a position angle of 226° covering the NE- 
Orion-S field, which passes close to the southeast corner of 
our field. They found densities no higher than 12000 cm~^, 
slightly lower than our determinations but consistent within 
the uncertainties. 

5.3 Chemical abundances 

In Fig. 7 we present the spatial distributions of ionic and to- 
tal oxygen abundances derived from CELs for the two fields. 
As expected, the and O^^ abundance distributions are 
rather similar to those of the [N ii]/Hq; and [O lll]/H^ line 
ratios shown in Fig. 5. However, unexpectedly, we find struc- 
tures in the spatial distributions of the total O abundance of 
both fields as we can see in Figs 7c) and 7f). The mean val- 
ues of the 0/H ratio at the BB and NE-Orion-S amount 
to 8.49=b0.03 dex and 8.48=b0.05 dex, respectively, which 
are in complete agreement with the typical oxygen abun- 
dances -8.50 dex- obtained by previous works in different 
parts of the Orion Nebula (Esteban et al. 1998, 2004; Bla- 
grave et al. 2006; Mesa-Delgado et al. 2009; Simon-Diaz & 
Stasihska 2011). In the BB field, the O abundance variations 
are within the typical error of about 0.10-0.15 dex, while 
in the NE-Orion-S the range is slightly wider, finding the 
lowest values at the southwest corner of the field, precisely 
where the O^^ abundances are lower and Te([0 III]) has a 
larger uncertainty, as we mentioned in §5.2. In Fig. 7, we 
can see that the areas with the highest total O abundance 
coincide in general with those showing higher ra- 
tio. Moreover, the structure of the spatial variations of the 
0/H ratio shown in Figs 7c) and 7f) is remarkably similar 
to that of the Ue maps of Figs 6a) and 6e). To understand 
this behaviour we have to consider that ionic abundance de- 
terminations derived from CELs depend very much on the 
physical conditions of the gas. 0^/H^ and 0^^/H^ ratios 
are both strongly dependent on Te, but not on Ue at the 
same level. In particular, the 0^/H^ ratio derived from the 
[O ll] 3727 A doublet is much more dependent on Ue than 
the O^^ abundance calculated from the nebular [O III] 4959, 
5007 A lines at least at the densities measured in our PMAS 
fields. Therefore, the striking similar structure of the 0/H 
and Ue maps may be produced by the strong dependence of 
the O^ abundance on Ue. The simpler explanation is that 
the ne([S ll]) we obtain for the ionization fronts is not indi- 
cating the true density of the O^ zone and, therefore, not 
appropriate for determining the O^ abundances. A possible 
way to check this solution could be to obtain high-spectral 
resolution spectra to measure and separate the [O ll] 3726, 
3729, 7319 and 7330 A lines and derive the precise physi- 
cal conditions of the O^ zone. Another explanation is that, 
even in the case that the physical conditions are the same 
for and O^, ne([S ll]) could not be the true one because 
the reported densities are about -specially in the case of 
NE-Orion-S- the critical densities of the [S ii] 6717, 6731 A 
lines -2000-3000 cm"^ and 10000 cm"^, respectively. In fact, 
considering densities somewhat higher than the true values 
produce higher 0^/H^ ratios. This is because we overes- 
timate the collisional de-excitation affecting the [O ll] line 



Table 2. Average distances^ from the peaks 
of the selected emission lines with respect to 
that of [O i] 





Bright Bar 


NE-Orion-S 


[S II] 


1(^0 ± 0'/4 




[O II] 


2'/3 ± 0'/9 


± I'll 


H I 


4(^5 ± V/0 


?>'.'! ± I'.'l 


He I 


5^5 ± (y/9 


3^3 ± 0^8 


[O III] 


r/6 ± v/1 


3'.' 9 di O'.'S 



" Each value is given as < > ±<t. 



fiuxes. Furthermore, [O ii] 3726, 3729 A emission fiux can be 
affected in the same way than [S ii] lines due to they have 
similar critical densities -3000 cm~^ and 5000 cm~^, respec- 
tively. Assuming that the main contribution comes from us- 
ing a wrong density, we have estimated by which amount we 
have to correct Ue in order to obtain the mean nominal value 
of the O abundance of the Orion Nebula -8.50 dex- in the 
zones of our fields with higher 0^/H^ ratio. This exercise 
gives that densities about 1000 and 4000 cm~^ lower than 
those determined for the BB and NE-Orion-S fields, respec- 
tively, are enough to wash out the problem. These values 
are rather modest and of the order or even smaller than 
the uncertainties estimated for ne([S ll]) and ne([Fe ill]) in 
each field (§4.1). Higher signal-to-noise ratio [Fe ill] spectra 
and/or good determinations of other indicators with higher 
critical densities such as C ii] 2326, 2328 A, [CI ill] 5517, 
5537 A or even C ill] 1907, 1909 A, as well as the observa- 
tion of [O ii] 7319, 7330 A to derive O^ abundances would 
be necessary to disentangle this situation. 

Finally, in Fig. 8, we present the O^^ and C^^ abun- 
dances derived from RLs for the NE-Orion-S, the only field 
where these emission lines were satisfactorily detected and 
measured. From the comparison of Fig. 7e) and Fig. 8a), we 
can see that the spatial distributions of 0^^/H^ determined 
from CELs and RLs show a rather similar behaviour, placing 
the maximum abundance values at the spaxels of the north- 
east corner, where the gas of higher excitation is located. 
Ionic abundance maps from RLs are nosier than those from 
CELs due to the faintness of the RLs. As it is common in 
ionized nebulae, we find that the O^^ abundances deter- 
mined from RLs are higher than those from CELs, which 
is related to the abundance discrepancy problem as we will 
discuss in §6.2 

6 DISCUSSION 

6.1 Structure of the ionization fronts 

In this section we study the ionization fronts of the two ob- 
served fields. With this aim, we analyse the spatial profiles 
of emission lines emitted from ions with different ionization 
potentials (0°, S+, 0+, H+, He+ and 0^+). We have ex- 
tracted the spatial profiles from cuts roughly parallel to the 
radial direction to ^^Ori C, connecting the southeast with 
the northwest corner in the BB field and the southwest with 
the northeast corner in the NE-Orion-S field. The use of all 
possible profiles along these directions allowed us to statisti- 
cally determine the relative distance between the maximum 
of the spatial distribution of the [O l] emission and the other 
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Figure 9. Spatial profiles of several emission lines traced along the main diagonal from the southeast to the northwest corner of the 
Bright Bar (top left) and along the third line below (to the southeast) of the main diagonal from the southwest to the northeast corner 
of the NE-Orion-S (top right). The zero point is located at the maximum of [O i] line emission and the offsets are given in arcseconds on 
the plane of the sky. Positive values in abscissae indicate decreasing distance with respect to ^-"^Ori C. The plotted lines are: [O l] (dotted 
line), [S l] (black solid line), [O ll] (dashed-dotted line), H I (thinner solid line). He I (three dotted-dashed line) and [O ill] (dashed line). 
The emission line fluxes are normalized to their maximum emission. In the middle and bottom panels, the spatial profiles of Te([N ll]) 
(solid line) and Te([0 ill]) (dashed line), as well as ne([S ll]) are presented. Te([N ll]) and ne([S ll]) are normalized to their respective 
maxima: 9480 K and 4760 cm-^ for the Bright Bar and 10390 K and 14200 cm-^ for the NE-Orion-S. 



emission lines across the fields. To estimate these distances, 
we have used the rebinned maps without applying any Gaus- 
sian smoothing filter as we did for the Ha images presented 
in Fig. 1. The rebinned maps do not change the relative po- 
sitions of the maxima and allow us to explore the changes 
in the profiles interpolating at subarcsecond scales. Further- 
more, these maps allow to use a total of 29 spatial profiles 
to calculate statistically significant averages. In Table 2 we 
present the average values of these distances and their stan- 
dard deviations. In Fig. 9 we present spatial profiles of the 
emission lines only for one representative cut: along the main 
diagonal in the case of the BB; and along the third diagonal 
line towards the southeast of the main diagonal in the case 
of the NE-Orion-S. In Fig. 9 we also show the temperature 
and density distributions along the same cuts. For all pro- 
files, we have located the zero point at the position of the 
maximum of [O l] emission, which is the farthest one from 
^^Ori C, as we saw in Fig. 4 and see now in Fig. 9. 

Theory about the ionization structure of nebulae can be 
found in Osterbrock & Ferland (2006). In our discussion be- 
low, we state that the stratification should be driven by the 
ionization energy of the atom. This is approximately true. 



but it should also be mentioned that the local ionizing radi- 
ation field in a nebula, which control the distribution of the 
ionic states of heavy elements, is determined by the opacity 
of the two most abundant elements, hydrogen and helium. 
In terms of the emission lines fluxes presented in Fig. 9, this 
implies that both [O i] and [S ii] arise near the zone where 
hydrogen becomes optically thick and the presence of these 
ions is combined with thermalized electrons coming from 
hydrogen photoionizations. [O ii] as well as [N ll] emissions 
arise from the zone where hydrogen is ionized and helium is 
still neutral. Finally, [O ill] emission comes from the near- 
est zone to ^^Ori C, where hydrogen remains ionized and 
helium is once ionized. In the spatial profiles presented in 
Fig. 9, the small differences with respect to that expected 
stratification are probably a combined effect of the geometry 
of each ionization front and the density dependence of the 
observed flux. 

The ionization structure of the BB has been well stud- 
ied in the literature {e.g. Hester 1991; Tielens et al. 1993; 
Walmsley et al. 2000; Pellegrini et al. 2009; van der Wiel 
et al. 2009). This bar is basically edge-on and its plane has 
an inclination of 7° with respect to the line of sight (Pel- 
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legrini et al. 2009), allowing to resolve its ionization struc- 
ture as we see in Fig. 9. The MIF is a thin layer located 
between the maximum emission peaks of [O l] and [S ll]. 
Pellegrini et al. (2009) indicate that the maxima of H2 and 
^^CO line emissions are detected at about 12^^ and 20^' be- 
hind the MIF position, respectively. As we see in Table 2, we 
found that the peak of the [S ll] emission is located about 
1" ahead [O l], which corresponds to a physical separation 
of 2.1x10"^ pc assuming a distance to the Orion Nebula 
of 436±20 pc (O'Dell & Henney 2008). The maximum of 
[O 11] emission, which also coincides with the maxima of 
[N II] (14.5 eV) and [Fe ill] (16.2 eV), is found about 2^3 
-4.9xl0~^ pc- with a standard deviation of 0^/9, implying 
that this zone is spatially less defined than that of [S 11] 
emission. The spatial profiles of the H i and He i emission 
lines are more extended. The mean distance of the peak of 
H I emission is 4'! 5 (9.5xl0~^ pc) with a standard deviation 
of 1" . In the case of He i emission, its maximum is located 
between the H i and [O ill] zones -as it is expected con- 
sidering the ionization potential of this ion (24.6 eV)- at 
5'! 5 (11.6 xlO~^ pc) with respect to the peak of [O i] emis- 
sion. The He l emission distribution represented in Fig. 9 
corresponds to the line at 6678 A, which permits to avoid 
problems related to self-absorption due to it comes from a 
singlet state. Finally, the [O III] zone reaches its maximum 
value at 7% (16.1x10"^ pc). 

In contrast to the BB, the ionization structure of NE- 
Orion-S has not been described in the literature. Taking into 
account the distances estimated in Table 2, the distribution 
of the spatial profiles shown in Fig. 9 and comparing with 
the results for the BB, it is clear that we can not resolve the 
ionization structure of the NE-Orion-S as well as of the BB 
because this front is much tilted with respect to the line of 
sight. In their figure 5b, Wen & O'Dell (1995) present the 
geometry of their Orion Nebula model along a cut with a 
position angle of 226°, which passes close to the southeast 
corner of our NE-Orion-S field. At a distance of 30^^ from 
the center of our field, the model by Wen & O'Dell (1995) 
indicates that the MIF has an inclination of about 60° with 
respect to the line of sight and after that the MIF remains 
rather parallel to the plane of the sky. The properties of the 
nebular structure to whom NE-Orion-S belongs, the Orion- 
S cloud, have been recently reviewed by O'Dell et al. (2009) 
showing that it should be located in the foreground with an 
unknown inclination as well as an uncertain placement in 
the three-dimensional space, but at a similar distance with 
respect to ^^Ori C as that of the MIF in the background 
(0.182 pc). The inclination angle of the NE-Orion-S have 
to be such to produce the observed distances between line 
emission maxima given in Table 2, where, on the one hand, 
the emission of [S ll] and [O ll] and, on the other hand, the 
emissions of H I, He I and [O III] are not well separated 
at our spatial resolution. We can compare the data of Ta- 
ble 2 for the BB and NE-Orion-S to obtain a first order 
estimation of the NE-Orion-S inclination. Assuming a pro- 
portional relation between the bar inclination with respect 
to the plane of the sky and the average distances with re- 
spect to the [O i] maximum emission of Table 2, we have 
estimated that the plane of the NE-Orion-S should be tilted 
about 48° ±13° with respect to the line of sight. This an- 
gle was calculated as a weighted average of the individual 
estimations for each emitting zone. In addition, it is impor- 



tant to note that the NE-Orion-S has a complex structure 
of shadowed regions and other features as we see in Fig. 1 
as well as in the new combined image of the Huygens re- 
gion obtained by J.R. Walsh using a drizzle technique to 
produce a mosaic with a resolution of O^.^l over the whole 
field^. Furthermore, the proximity and unknown exact po- 
sition with respect to the Trapezium cluster is also playing 
a role on the geometry of what we observe as we mention 
in §5.1. Therefore, it can be difficult to interpret this ion- 
ization front as a simple slab with a given inclination. In 
this sense, our analysis of the individual spatial profiles ex- 
tracted parallel to the main diagonal of the NE-Orion-S field 
can give further information about the complex structure of 
the NE-Orion-S. We have noted that the emission of [S 11] 
and [O 11] are much separated in the profiles extracted be- 
low -to the southeast of- the main diagonal than above it. 
In fact, we have estimated that the average distances be- 
tween the emission peaks of [S 11] and [O ll] and that of [O i] 
are of about 0V7±0V4 and 1^.^9=b0^.^9 below the main diago- 
nal, and 0'.'1±0'.'3 and 0'.'6±0'.'3 above it. This fact suggests 
that the inclination of the NE-Orion-S is not constant along 
the observed field. Below the main diagonal, the plane of 
the NE-Orion-S seems to be less inclined with respect to 
the line of sight -more edge-on. But, above the main diago- 
nal, the inclination angle is larger and the peaks of the line 
emission of the different ionic species can not be resolved. 
The variation of the inclination of the NE-Orion-S along the 
field may also explain the change of the width of the area of 
high [N 11] /Ha ratio -a good indicator of the inclination of 
an ionization front, as it is commented in §5.1- across the 
field (see Fig. 5), which seems to be wider below the main 
diagonal. 

In Fig. 9 we have also plotted the spatial profiles of 
physical conditions of the BB and NE-Orion-S along the 
same diagonals as in the case of upper panels. Here, some 
of the results commented in §5.2 are more clearly seen: 
Te([N 11]) is always higher than Te([0 ill]), but they tend to 
be more similar inside te bars. For the NE-Orion-S we can 
note a correlation between the Ue and Te([N ll]) profiles, 
where lower temperatures are reached where the densities 
are higher. At the BB, this effect is not as clear as in the 
case of NE-Orion-S, this can be due to the lower densities 
of the BB. 

6.2 Abundance discrepancy problem 

The Orion Nebula is an unique object to analyse and explore 
particular problems in common with other H II regions. This 
is the case of the well known abundance discrepancy (AD) 
problem, i.e. the disagreement between the abundances of 
the same heavy-element ion derived from CELs and RLs. 
The AD problem is quantified making use of the AD factor, 
ADF, which can be defined for a given ion as 

ADF(X-) = /o,(|;)^^-/o,(|;)_. (2) 

From intermediate and high resolution spectroscopy, 
several works have found that the ADF(0^^) is always be- 
tween 0.1 and 0.3 dex in Galactic and extragalactic H 11 

^ J.R. Walsh webpage: http://www.stecf.org/ jwalsh/ 
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Figure 10. Spatial distribution of ADF(02+) for the NE-Orion-S 
with Ho; contours overplotted. 



regions (see e.g. Esteban et al. 2002; Tsamis et al. 2003; 
Esteban et al. 2004; Garcia-Rojas & Esteban 2007; Lopez- 
Sanchez et al. 2007; Mesa-Delgado et al. 2008; Esteban et al. 
2009; Mesa-Delgado & Esteban 2010). Garcia-Rojas & Es- 
teban (2007) argued that the AD problem in H ii regions 
seems to be consistent with the predictions of the tempera- 
ture fluctuation paradigm proposed by Peimbert (1967) and 
parametrized by the mean square of the spatial variations 
of temperature, the so-called parameter. Under this hy- 
pothesis, the AD would be produced by the very different 
temperature dependence of the emissivities of both kinds of 
lines and the abundances from RLs would be the represen- 
tative ones of the ionized gas. However, the existence and 
origin of such temperature fluctuations are still controver- 
sial and the AD problem remains open. Other authors as 
Tsamis & Pequignot (2005) and Stasihska et al. (2007) have 
proposed a different hypothesis to explain the AD problem 
in H II regions. Based on the model for heavy-element mix- 
ing of Tenorio-Tagle (1996), Stasihska et al. (2007) propose 
that the presence of cold metal-rich droplets of supernova 
ejecta still not mixed with the ambient gas of the H ii regions 
would produce most of the RL emission, while GEL emission 
would be produced by the ambient gas, which would have 
the expected composition of the local interstellar medium. 
In such case, the abundances from GELs would be the truly 
representative ones of the nebula. Recently, the studies per- 
formed in the Orion Nebula by Mesa-Delgado et al. (2009), 
Tsamis et al. (2011) and Simon-Diaz & Stasihska (2011) 
have provided new observational clues to address this puz- 
zling issue. Mesa-Delgado et al. (2009) obtained echelle spec- 
troscopy of the HH 202-S separating the kinematic com- 
ponent of the background emission from that of the gas 
flow finding clearly higher ADF(0^^) values at the high- 
velocity component. This result suggests a possible relation 
between the ADF(0^^) and the gas flow velocity. Tsamis 
et al. (2011) carried out a high-spectral resolution analysis 
of the protoplanetary disc LV2 with FLAMES. Their study 
of LV2 emission with the nebular background subtracted 
shows ADF(0^^) values which are basically zero. This sur- 
prising result seems to be related to the high electron densi- 
ties expected in these kind of objects, which affect the ionic 
abundances derived from GELs. Therefore, it is necessary to 
explore the behaviour of the ADF(0^^) in other morpholog- 



ical structures of the Orion Nebula -as stratified ionization 
fronts- in order to better understand and constrain the AD 
problem. Finally, we consider important to comment the re- 
cent remarkable result by Simon-Diaz & Stasihska (2011), 
who find that the nebular oxygen total abundances derived 
from RLs agrees better with the stellar abundances of B-type 
stars of the Orion Association than those determined from 
GELs. This result is an indication that nebular abundances 
obtained from RLs may be more reliable than those derived 
from the standard method based on GELs and therefore 
stressing the importance of investigating the AD problem. 

As we pointed out in §5.3 we derived O^^ abundances 
from RLs in the NE-Orion-S field. Subtracting these abun- 
dances from those obtained from GELs (Fig. 7) we find the 
spatial distribution of the ADF(0^^), which is presented in 
Fig. 10. As we can see, the ADF(0^^) seems to be slightly 
higher at the northeast corner of the field, but the map 
does not show any remarkable structure inside the ioniza- 
tion front. The range of the variations goes from 0.05 to 
0.25 dex with a typical error of about 0.12 dex for each 
spaxel. However, we have checked that the extreme values 
of the ADF(0^^) correspond typically with spaxels with low 
signal-to-noise ratio in the O ii RLs of their spectra. There- 
fore, not particular variations of the ADF(0^^) are related 
to the presence of stratified ionization fronts, at least with 
the properties and physical conditions of the NE-Orion-S 
and the spatial resolution of our observations. The mean 
value of the ADF(0^^) in the field is about 0.14 dex with a 
standard deviation of about 0.05 dex. Therefore, there is a 
real difference between the O^^ abundances from RLs and 
GELs, which is in the range between 21% and 44%. This 
is precisely the value found in our previous studies of the 
Orion Nebula {e.g. Esteban et al. 2004; Mesa-Delgado et al. 
2008). 



7 CONCLUSIONS 

In this paper we analyse integral field optical spectroscopy 
of two two remarkable structures related to ionization fronts 
of the Orion Nebula: the Bright Bar (BB) and the northeast 
edge of the Orion South cloud (NE-Orion-S). In both fields, 
we have mapped emission line fluxes and ratios, physical 
conditions and chemical abundances at spatial scales of l^^ 
The maps clearly show the ionization stratification of 
both ionization fronts, which present rather different incli- 
nation angles as well as different distances with respect to 
^^Ori G, the main exciting source of the nebula. The ioniza- 
tion structure of the BB is well resolved because the ioniza- 
tion front is almost edge-on -inclination angle of 7° with re- 
spect to the line of sight-, while the NE-Orion-S is a portion 
of the Orion-S cloud suspended within the main body of the 
nebula, which is more tilted and presents a more complex 
structure. We have analysed spatial profiles of selected emis- 
sion lines of ions with different ionization potential and esti- 
mated the average distances of the peaks of the emission of 
these lines with respect to the maximum of [O i] . Gomparing 
those distances in both fronts, we have estimated that the 
plane of the NE-Orion-S field could be tilted about 48° ±13° 
with respect to the line of sight. However, the distances be- 
tween the peaks of the lines with respect to [O i] change in 
the different spatial profiles extracted for the NE-Orion-S, 
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indicating that the incUnation angle of this ionization front 
should vary across the observed field. 

The maps of electron density, ne(S ll]), show that NE- 
Orion-S is much denser than the BB. The NE-Orion-S has a 
mean density of about 10900 cm~^, while the BB has 4370 
cm~^. In both cases, the [S ll] lines may be somewhat af- 
fected by collisional de-excitation because the nominal val- 
ues of Ue we find are of the order of their critical densities. 
In any case, this de-excitation should not be severe due to 
the rather similar Ue we determine from [Fe ill] lines. The 
structure of the temperature maps are rather featureless and 
Te([N ll]) is usually higher than Te([0 III]), specially outside 
the bars. The maps of total O abundance show mean values 
consistent with previous determinations of the literature for 
the Orion Nebula, but they show a pattern very similar to 
the maps of 0^/H+ ratio and ne([S ll]). The spaxels with 
higher Ue values also show higher and O abundances. We 
find that considering values of Ue slightly lower than those 
determined from the [S il] line ratio -but of the order of 
the quoted uncertainties- the problem can be satisfactorily 
solved. 

Finally, the higher signal-to-noise ratio of the spectra 
in the NE-Orion-S field allows us to detect, measure and 
map the faint RLs of O ii and C ii. We explore the abun- 
dance discrepancy problem in this field comparing the O^^ 
abundance maps derived from CELs and RLs, finding a fea- 
tureless spatial distribution for the abundance discrepancy 
factor of O^^, ADF(0^^), which shows values ranging be- 
tween 0.05 and 0.25 dex with a typical error of about 0.12 
dex. The mean ADF(0^+) is 0.14zb0.05 dex, in excellent 
agreement with other previous determinations in other areas 
of the Orion Nebula -except at Herbig-Haro objects (Mesa- 
Delgado et al. 2008, 2009) and the LV2 protoplanetary disc 
(Tsamis et al. 2011). This result implies that the O^^ abun- 
dances from RLs are between 21% and 44% larger than those 
from CELs and that there is not any particular variation of 
the ADF(0^^) related to the presence of stratified ionization 
fronts, at least for the properties and physical conditions of 
the NE-Orion-S and the spatial resolution of our observa- 
tions 
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